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Abstract 

An application of approximate transformation groups to study dynamics of a sys- 
tem with distinct time scales is discussed. The utilization of the Krylov-Bogoliubov- 
Mitropolsky method of averaging to find solutions of the Lie equations is considered. 
Physical illustrations from the plasma kinetic theory demonstrate the potentialities 
of the suggested approach. Several examples of invariant solutions for the system 
of the Vlasov-Maxwell equations for the two-component (electron-ion) plasma are 
presented. 

1 Introduction 

In analyzing various physical systems we frequently deal with a situation when a com- 
plicated dynamics of different systems appears as a superposition of "fast" and "slow" 
motions with incommensurable characteristic scales, for example, slow evolution of "back- 
ground" system characteristics accompanied by fast oscillations in the vicinity of a back- 
ground state. This type of behavior seems typical for various linear and nonlinear prob- 
lems (numerous examples are found in [TJ[2]), e -S- f° r celestial mechanics in studies of a 
motion of planets, for mechanics when treating oscillatory regimes of systems with slowly 
varying parameters, for various nonlinear problems of mult i- component plasma. 

The availability of different scales (though the origin of these scales depends upon the 
particular system of interest) allows to simplify the analysis of the complicated dynamics 
by treating "fast" and "slow" motions separately. These ideas underlie an essence of 
asymptotic analytical approaches, the method of averaging IHlH], the method of multiple 
scales pp, and other asymptotic methods (see, e.g., [HE])- 

As to relation of modern group analysis to nonlinear dynamics here we will point 
to an interpenetration of ideas from both fields: on one hand the use of the Lie group 
theory in asymptotic methods for integration of nonlinear differential equations gives (in 
combination with the Hausdorff formula) the theoretical basis for the method of averaging 
[3] and provides a regular procedure for calculating the asymptotic series in this method 
jSJElE]. On the other hand an introduction of multiple-scales approach to modern group 
analysis [7] enhances the potentiality of approximate transformation groups [8]. 
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The present paper uses the Krylov-Bogoliubov-Mitropolsky method (KBM-method) 
of averaging in group analysis of the system of equations that describes the evolution 
of plasma particles in multi-component plasma. The linearity of the group determining 
equations plays the decisive role in separating fast and slow terms in coordinates of a 
group generator and in successive use of the KBM-method for constructing the asymptotic 
solutions of the Lie group equations. 

The paper is organized as follows: in Section [2] we introduce the basic equations for 
evolution of a two component plasma and introduce small parameters which give rise to 
different time and spatial scales for distinct plasma components. In the next Section [3] 
we describe the solution of determining equations, which define approximate Lie point 
symmetry group. The use of KBM-method in finding solutions of the Lie group equations 
constitutes the basis of the Section HI Several examples of an application of the suggested 
approach are presented in Section \5\ and Sectional In the Conclusion we discuss the results 
obtained and the future application of the KBM-method in modern group analysis. 



2 Basic Equations: Electron-ion Plasma 

We start with kinetic equations for distribution functions, f e and /*, 

d t f e + v e d x f e - (e/m)Ed v ef e = , d t f + v%f + {Ze/M)Ed v ,f l = , (1) 

for both species of two-component plasma consisting of electrons and ions with mass m 
and M and charges e e = — e and e l = Ze, where Z is a charge number, and equations for 
a self-consistent electric field E 



d x E = Airp , d t E = -Attj , d t p + d xJ = . 



(2) 



Here charge p and current j densities are related to moments of the distribution functions 
via nonlocal material relations: 



P 



dv l f 



dv e f 



J 



dvVf 



dv e v e f 



(3) 



Equations (DO)-© are known as a system of the Vlasov-Maxwell equations for a collisionless 
plasma. We are interested in the solution of the Cauchy problem for kinetic equations ([T]) 
with the initial conditions 



t=o 



fo(x,v e 



r 



t=0 



(4) 



which depend on a particular physical problem. In what follows we consider an evolution 
of localized plasma bunches and assume sufficiently smooth (e.g., Maxwellian) initial 
distribution functions with electron T e and ion Tj temperatures and initial densities of 
electrons n e (x) = j dv e fy and ions = j dv l f$ with the characteristic scale L. Below 
we consider a typical situation when L is much greater than the Debye radius of electrons 
r De = \jT e / {Annie 2 ). The difference in mass of plasma particles specifies two different 
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time scales, namely dimensionless time uj^t for "fast" electron motions and r — fit for 
"slow" motions, fi = a/ Zm/M C 1. It is natural that electrons are involved in both 
fast and slow motions, hence the electron distribution function depends on both t and 
r. On the contrary, we consider ions not involved in fast motions. It means that the 
ion distribution function does not depend upon fast time t, but only upon slow time 
r. It also means that averaging upon fast time eliminates the fast component E of 
the electric field E = E + E in the ion kinetic equation that contains only the slow 
electric field E. Then introducing dimensionless variables, electron velocity u = v e /Vr e , 
Vr e = yTefrn, dimensionless ion velocity w = v l /c s , c s = a/ ZT e /M, dimensionless 
electric field p = e(eEL/T e ) , e = r De /L <C 1 and dimensionless distribution functions 
f e = ( n o/^ / Te)g, P = ( n o/(Zc s ))f , = n e W(0) we come to the following system of 
basic equations in dimensionless variables 

d t g + /j,d T g + eud x g - pd u g = , d T f + ewd x f + pd w f = , (5) 

ed x p = I dwf — / dug , d t p + \id T p = —ji / dwwf + / duug . (6) 



From the group analysis point of view this system of equations should be supplemented 
by the four additional equalities, 

d w g = , d u f = , d u p = , d w p = , (7) 

which are evident from the physical point of view. 



3 Lie symmetry group 

The Lie point symmetry group admitted by the system fl5]) and fl6]) is defined by a sym- 
metry group generator 

x = eo t + ed x + ed u + + eo T + v% + + v 3 d f . (8) 

In the canonical form this generator is written as 

Y = K x d g + K 2 d p + K 3 d f , (9) 

k 1 = v 1 - eo t g - eo x g - eo u g - eo T g , 

k 2 = rf - ^d t p - i 2 d x p - £ 5 d T p , 

= v 3 - eoj - ed w f - ed T f . 

When applying group generator ([9]) to ([5]), and ([7]) we get a system of determining 
equations for coordinates £\ rf of the generator ([8]), 

DtK 1 + [iDt-k 1 + euDxK 1 — pDuK 1 — n 2 d u g = , 

D t k 3 + ewD x K 3 + pD w K 3 - R 2 d w f = , (10) 
D w k x = , D u k 3 = , D u k 2 = , D w k 2 = , 
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eD x k 2 - / dwK 3 + / dun 1 = , 

> f f (U) 

D t n 2 + fxD T K 2 + fi / dwwn 3 — / duuK 1 = , 

which should be solved in view of the basic equations ([5]) , (J6]) , (JTJ) and all their differential 
consequences. Here D t , D T , D x , D u and D w are total differentiations with respect to the 
variable denoted by lower index, 

D t = d t + (d t g)d g + (d t p)d p , 
D T = d T + (d T g) d g + (d T f) d f + (d T p) d p , 
D x = d x + (d x g) d g + (d x f) d f + (d x p) d p , 
D u = d u + (d u g) d g , D w = d w + (d w f) d f . 

To find coordinates £ J , rf from the system of determining equations we use the approach 
developed in [9] (see also [T5l Chapter 4]). Following this technique we separate the 
determining equations for and rf into local determining equations, (ITOl) . which arise 
from invariance of ([5]), (JTj), and nonlocal determining equations, ( Till , which follow from 
invariance conditions for ([6]). Solutions of local determining equations give the so-called 
intermediate symmetry [9]. 

Two distinct moments should be taken into account here: first, in view of multi-scale 
dynamics we outline in coordinates £ 2 , £ 3 and rf, entering local determining equations 
the fast terms denoted by variables with tilde and slow terms denoted by variables with 
bar 

C-C-C; 2 = 1,2,3, r j 2 =f + f. (13) 

Due to the fact that both local and nonlocal determining equations are linear in £ and rj, 
we thus can separate terms of different characteristic scales. Then omitting trivial tedious 
computations we rewrite fast 

rf + V (dtl 1 + /^rf 1 - d u ^) + d t i 3 + jid T ~e + eupdj 1 + eud x £ 3 = , 

. \ _ / ; .x _ (14) 

e^ 3 - dti 2 - jidri 2 + eu (dti 1 + jidrZ 1 - d x i 2 ) + e 2 u 2 d x ^ = , 

and slow local determining equations, 

fj 2 +p (<9 T £ 5 + ewd x £> - d w ^) - d T t - ewd x ^ = , 
rf + p (jid T ? - d u £ 3 ) + jid T ? + eupd x ? + eud x ? = , 

- l&r? + eu (jidrl 1 - d x l 2 ) + e 2 u 2 d x l l = , 

- d4 2 + ew (d T e - d x ?) + e 2 w 2 d x e = . 

Here, in f[T4"l) and (fT5|) the dependencies of £*, £\ and rf, rf upon group variables are given 
by 

f =i\t,T,x), £ 2 = i 2 (t,T,x), £ 3 = £ 3 (t,T,X,u), f = fj 2 (t,T,X,p), 

? = t\r,x), ? = ?(t,x), C ! ?(t,x,u), ? = t\r,x,u), (16) 
£ 5 = /^\ V 2 =f(T,x,p), r, 1 = r, 1 (g), T 1 3 = rf(f). 
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Second, we shall take an advantage of small parameters in (jSJ), ([6]) and, as is customary 
in the approximate group analysis technique [8], express the coordinates of the group 
generator as power series in e and /x, 



oc 



C=J2 eVf (W) , rf = eW ( *'° ■ (17) 

k,l=0 k,l=0 

Collecting terms of the same order we come to the following infinite set of equations that 
relate coordinates of different orders for the fast 

v mi) + v (d£ 1(k ' l) + (i - Wi^-v - fi^fW)) 

+ d t i 3 ^ + (1 - 5 h0 )d T i 3 ^ + up(l - S^dJ 1 ^ 
+ u(l - 6 kfi )d x i 3 ^-^ = , k,l^0, 

-d t ~e [Kl) + (i - W^*" 1,0 - (i - 

+ «(1 - 4,o) (ftf^ + (1 - 5 lfi )d T ^ k -^ - fi^f (*-^ 
+ (i-4,o)(i-4,iy^| 1(fe - 2 ^ = o, 

and the slow terms 

ffiH) + p (^(V) + (1 _ 4 i0 ) w ^( fe -M) _ ^4(M)) 

- ^ 4(M) - (i - 5 kfi )wd x e {k ~ l > l) = o , 

V 2{k ' l) + P {(1 - 5 l>Q )d T ^ k > l -V - fi^»(*.0) + (i _ 5 z , )<9 r £ 3 ^ 
+ (1 - ^up^*- 1 ^ + (1 - ^,o)«^f (fe - 1 ' Z) = , 

(i - 4,o)f (fc " 1 '° + (i - 5 kt0 )u ((i - ^o)^*- 1 - 1 - 1 ) - d x ^ k -^) 

- (i - M^C 2 ^ + (i - 4,o)(i - 4,i)« 2 ^e 1( ^ 2,0 = , 

(1 " 5 fe) o)e 4(fc - M) - Or?™ + (1 - (flr^*- 1,0 - d*¥ {k ' l) ) 

+ (l-5 kt0 )(l-5 kil )w 2 d x e {k ~ 2 ' l) = 0- 

Using approximate intermediate symmetry, which follows from solutions of equations 
(}T5|) - (jl9p . in non-local determining equations (ITTj) we find a solution of the latter using 
variational differentiation (see j9], [151 ch.4] for details) and obtain the symmetry of the 
complete system (jSJ), © and ([7j). 

In constructing the solution of the b.v.p. (SJ) — (ED we need not the entire set of genera- 
tors but rather such a combinations of group generators that leaves invariant the pertur- 
bation theory solution in powers of t and r, the so-called renormalization group symme- 
tries [10J. Hence, we should specify the initial particle distribution functions f° = f\t=o, 
g° = g\t=o- For concreteness, we assume the initial velocity distribution functions to be 
maxwellian: 

g° = n e (x) exp(- M 2 /2) , f° = n» (ar) exp(- W 2 /2 7 2 ) , (20) 



(19) 



5 



with the initial densities tIq(x) and n^x) and the initial zero average velocities. In account 
of these initial distribution functions we have the following initial electric field 

p°(x) = (1/e) / dx (ni(x) - n e (x)) . (21) 
Jo 

Perturbation expansion of the Cauchy problem solutions in powers of t and r gives terms oc 
0(t) and oc 0(t) for the electron distribution function and oc 0(t) for the ion distribution 
function, and oc 0(t 2 ) and oc 0(t 2 ) for the electric field. Invariance conditions for these 
solutions specify the coordinates (fT7|) of the group generator ([8]). Leaving only terms that 
are linear in e and p we write these coordinates as follows 

f 1 = 1 + E T 2 d x £ , £ 2 = e ((5/n) sin fit + //t£) , 

£ 3 = 5 cos fit — Eprud x ^ , £ 4 = /i (£ — erwd x £) , 

r? 2 = 5fi sin fit - Sperpd^ , ^ 5 = /i^ 1 , (22) 

^=^ = 0, n 2 =n i (r,x) = J dwf . 

The dependence of functions £(x) and S(x) upon x is expressed in terms of the initial 
densities distributions and the initial electric field p°, 

£ = -e ({d x nl/nl) + 7 2 (^n l /n l )) , 6 = -p° - e(d x n e /n e ) , 7 = V Ti /c s . (23) 

For arbitrary parameters e and p, and arbitrary initial density distributions n[j' % formulas 
(|22|) describe the approximate symmetry. However, in two limiting cases infinite series 
(fTTl) terminate and we get the exact symmetry group. The first case is referred to electron 
plasma with neutralizing homogeneous ion background (p = 0, fi 2 = n = const) [TTHTB"] . 
which gives the generator 

X = d t + e (5 /SI) sin Vttd x + S cos Vttd u + 80, sin fit<9 p . (24) 

The second case is referred to quasi-neutral approximation for electron-ion plasma with 
zero current and charge densities j = p = [2] that is realized for 5 = and the initial 
gaussian densities distribution, £ a (3x 

X = (l + (3et 2 ) d T + ep>Txd x + (5 (px — etu) d u + (3 (x — etw) d w — 3ef3rpd p . (25) 

The additional term in £ 3 in (125]) that refers to acceleration of electrons is omitted in (1221) 
as it is of the higher order 0(p 2 ) as compared to that included in (1221) . 
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4 Lie equations 



To construct group invariant solution for the b.v.p. (SJ) — (]6j) we should find solutions of 
the Lie equations for the group generator (jSJ) with coordinates ([Tj 



dt 

da 
dx 

da 
du 

da 
dw 

da 
dr 

da 
dp 

da 

df dg , 

-1- — -J- — U ) 9\a=0 = g , J\a=0 = J ■ 

da da 

Solution of the b.v.p. (J4j) — (JHJ) are expressed as usual in terms of invariants of the group 
dHJ), ( 1T71) that result from solutions of ( 1261) after excluding the group parameter a. Due to a 
difference in characteristic time scales we can separate "fast" and "slow" group invariants, 
applying the averaging procedure to Lie equations. 

In fact, at small time t > 0, ^> t ^> I/O, the "ion" terms that are oc fi can be 
omitted and we come to simplified Lie equations (equations for group invariants /, g, r, 
w are omitted here), 



1 + st d x £ , t\ a=0 = t' , 

e {{S/fl) sin Qt + /ir^) , x\ a=0 = x' , 

5 COS Qt — EfiTudxC, , M| a =0 = u ' j 

= fi (£ - £Twd x £) , w\ a=0 = w' , (26) 
H (l + £r 2 9 z ^) , T| a=0 = r' , 
50 sin Ot - ZijLETpd^ , p| a=0 = p' , 



ctb dx 

— = 1, t| a=0 = t'; — =e{S/Q)smQt, x\ a=0 = x' ; 
da da 

— = 5cosQt, u\ a=0 = u ; — = 5QsmQt, p\ a =o — V , 
da da 

the solutions of which define invariants of "fast" motions at small time £ <C 

■h = P + ScosQt = -e(d x nl/nl)\ x=x i , 

J 2 = x + (e5/n 2 ) cosOi = x' + (e5(x')/n 2 (x')) , 

J 3 = u — (5/0) sinOt = u , 

Ja = g = g°(x',u) . 



(27) 



(2f 



On the contrary, averaging the complete Lie equations on a large time scale T > we 
come to Lie equations defining "slow" motions (equations for group invariants are again 
omitted), 

dr 2 - rfx - dp - 

— = 1 + er 9^ , — = er£ , — = -3erpd s £ , 



with the corresponding "slow" invariants 



h = f = f, l2 = g = g', h=P?=p't' 3 



2 fx /2 



Here the "primed" variables are related to values at r — )■ 

/' = / (sX), g^g'M, p' = -ed sV , s = x'+ £ -§fl 

U z {x') (31) 
e 2 ^^ + 7io(ar) - e v = , 9 b ^|«b=o = §»¥>|a>-+oo = , p\ x=0 = C < oo. 

In the next section we use the fast and slow invariants to construct analytical solutions 
of the Cauchy problem for the kinetic equations ([ll -p]) . 

5 Slow dynamics of plasma particles 

Let we consider the slow dynamics of plasma particles under simplifying assumptions, 
small value of p < 1, and low ion temperature, 7 — > 0. Then, following (13"0 ]) - (f3"T]) . s~x' 
and £ coincides with p' = p\ T= o, and we come to simplified expressions, which define 
dynamics of plasma ions, 

f = f°(xW), g = g°(x',w'), p=(£'/0 3 pV), st 2 = 2? [ X dz/e, 



x' 



1 



y 



I dy( I dz/e) , 



-1/2 



(32) 



For completeness we also present global characteristics for plasma ions, their average 
velocity v l av , density n l av and temperature T\ 

*• = m £ d « {£ v) ^ • "™ = • T ' = Ti (f ) 2 • (33) 

These formulas are analyzed below for two distinct initial electric field and density distri- 
butions. 

5.1 Examples of slow plasma dynamics: gaussian density profile 

We start with a specific situation when electron and ion densities balance each other, and 
are described by gaussian curves, 

n° e (x) = n°(x) = (1/v 7 ^) exp(-a; 2 ) , £ = lex (l + 7 2 ) . (34) 
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Substituting expressions ( pUD in fl32|) - fl33|) we get formulas for finite group transformations 
for "slow" variables 

P= n ~ P 2 2W2 , x = xWl + v 2 r 2 , z/ 2 = 2e 2 (l + 7 2 ) , 

( + } , / 2 x (35) 



u = = , w 



and formulas for ion average density, velocity and temperature 



w'-\ x't , f = f , g = g' . 



1 fx 



^u{1 + v 2 t 2 ) V 1 + 



T.2 



= ^ZT^fTfV eX P I i.l./W, 

(36) 



x t/ 2 r „, T(5 



^ £ (1 + v 2 r 2 ) ' T (1 + ^2)* 

These formulas demonstrate the self-similar dependence of the ion density - the Gaussian- 
type density distribution is preserved though the spatial scale of this distribution as well as 
the maximum of the ion density varies. The spatial dependence of the average ion velocity 
is linear, and the ion temperature is uniform in space and monotonically decreases with 
the growth of time r, as it was demonstrated in [TJ]. In case, when the initial ion density 
varies from the electron density, i.e. for n®(x) = (b/ y/n)e~ b2x2 , b ~ 1, the above formulas 
are still valid, provided 1 + 7 2 is replaced by 1 + 6 2 7 2 . The difference between the initial 
densities distributions of particles leads to nonzero values of p°, 

p% = 1 (Erf (fox) - Erf (x)) , £ ga = 2ex (l + fo 2 7 2 ) , 5° = - P ° gs + 2ex . (37) 

however for b close to unity the values of p° are small as compared to £. Figure [T] (left 
panel) illustrates this fact, while the right panel shows the spatial distribution of the 
"slow" electric field p at different nonzero time moments. 

As for the spatial dependencies of the average ion velocity and density at the same time 
moments they are presented on Figure El It follows from ([37j) that the oscillating electric 
field is of the order of the average electric field that accelerates ions. In the next section 
we consider the opposite situation when the initial electric field practically concise with 
the initial "slow" electric field p and the amplitude of the "fast" electric field is small. 

5.2 Examples of slow plasma dynamics: Lorentz density profile 

Turn now to the case when there is very slight difference between p° and £ that is realized, 
for example, for the Lorentz-type initial density distribution, 

n° e (x) = (1/tt(1 + x 2 )) , n?(x) = (6/tt(1 + foV)) , |fo-l|<l. (38) 

Substitution of (138]) into (T2Tj) gives the following formulas for the spatial distribution of 
the initial electric field and the function £, 

1 (I 7 2 fo 2 \ 

p?(x) = — (arctan for — arctan x) , £i = 2ex I H — — ) . (39) 

ne \ 1 + x 2 1 + b 2 x 2 J 
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Figure 1: Distribution of electric field p° gs (blue line) and £ gs (red line) for r =0 (left), and 
"slow" electric field p (right) for r = 2(blue line), r = 8(red line), r = 12 (yellow line), 
t = 18(green line), and for b = 1.01, e = 0.1, fi = ^1/2000 and 7 = 0.1. 




Figure 2: Averaged density and velocity distributions of ions for r = (light-blue line), 
r = 2(dark-blue line), r = 8(red line), r = 12 (yellow line), r = 18(green line), and 
b = 1.01, e = 0.1, fi = v/1/2000, and 7 = 0.1. 



The left panel of Figure [3] demonstrates the difference between p° and £j, and the right 
panel shows the spatial distribution of the "slow" electric field p at different moments of 
time t. 

As for the average ion density, temperature and velocity they are given by the formulas 
h C ' Z* = 2s(tY\ < = Az= I* dy ( [ V dz/tfY 1 ' 2 . (40) 



av 7r(i + 6V 2 )6' °V6; ' ^ ZiV* 

and are plotted on the Figure HI 



6 Fast dynamics of particles 

In this section we use slow invariants to restore the complete dynamics of fast particles, 
electrons. For clarity's sake we consider the case of small values of 5e 1, which means 
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Figure 3: Plots of p° (blue line) and £z (red line) at r = (left), and "slow" electric field p 
(right) at r = 4 (blue line), r = 10 (yellow line), r = 18 (red line), for a = 1.0, b = 1.0661, 
£ = o.l, ^ = v /I/2000 and 7 = 0.001. 





Figure 4: Density and velocity distributions for ions at r = (green line), r = 4 (blue 
line), r = 10 (red line), r = 18 (yellow line), for b = 1.0661, e = 0.1, fi = ^1/2000 and 
7 = 0.001. 

that x is identical to x, and rewrite the Lie equations f|26|) in a simplified form 

^ = l + er 2 d x ^ ^2 = ^ cos fit, ^M = ^ 3 cosfit. (41) 
da da da 

According to the procedure of averaging [3] we can write the solutions of equations (14"T!) 
by integrating over fast time t and taking into account the dependence upon slow time by 
including the dependence upon x and r into 5 and £. However, the enhanced precision is 
achieved by direct integration of the Lie equations (j4Tj) in account of the slow dependence 
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of r upon x as given by slow motion invariants 



P £3 P + £3 J d °° efJLT ( x 



u 



Cw' + i / dx" 



(42) 



Efxr(x 



-co S (Qt(x")). 



Electron distribution function g = g' \ 
Thus substituting x = x and u' from 
the integral characteristic, the average electron velocity and density 



g°(x', u') is the invariant of group transformations. 
^2]) in g° and integrating over the velocity u gives 



n ;, 



n e n (x 



') (eve) 



i 



dx" 



e\it{x" 



cos(m(x")) 



(43) 



To illustrate these formulas we employ results of the previous section and consider the 
Lorentz-type initial densities profiles (j3~8j) with the function £ = £z defined by ( |39l) . Sub- 
stituting £/ in (l42l)-( 143|) we obtain the formulas for the spatial distribution of the electric 
field and the average electron velocity that are presented on the figures below for three 
different time moments. Figure [5] corresponds to moderate values of r = 4, when the 
ion density is concentrated mainly in the center of the bunch thus leading to small-scale 
spatial oscillations primarily in this region. As the bunch spreads with growth of r the 





Figure 5: Electric field and average electron velocity distributions at r 
b = 1.0661, e = 0.1, n = ^1/2000 and 7 = 0.001. 



4 for a = 1.0, 



small-scale spatial oscillations moves outward as shown on Figured for r = 10, and on 
Figure [7] for r = 18. The figures also show that the "complete" electric field p in this case 
oscillates with the same spatial period as the mean electron velocity and only slightly 
differs from the average electric field p (compare with Figures [3]). 
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Figure 7: Electric field and average electron velocity distributions at r = 18 for a = 1.0, 
b = 1.0661, e = 0.1, fi = v/1/2000 and 7 = 0.001. 

7 Conclusion 

In the above analysis of a particular physical problem, expansion of a plasma bunch, a 
new promising tool for analyzing nonlinear multi-scale systems was considered. The main 
idea consists in employing the Krylov-Bogoliubov-Mitropolskii procedure of averaging 
to construct solutions of Lie equations. The procedure of separating "fast" and "slow" 
terms in coordinates of group generator naturally occurs in determining equations while 
constructing the symmetry for nonlinear equations that describe multi-scale behavior of 
any physical system. In our consideration we use the averaging procedure in combination 
with a perturbation technique of group analysis [8J that gives approximate symmetries 
for the analyzed problem and helps to construct approximate RG-invariant solution for 
arbitrary initial distribution functions of particles. 

The use of the averaging procedure in modern group analysis naturally separates in- 
variant manifolds, related to slow and fast Lie equations into slow and fast invariant man- 
ifolds. This separation is in the root of the theorem of invariant representation [16, §18]: 
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averaging the fast invariant solution that appears as an oscillating "curve" on fast mani- 
fold yields a smooth curve on the slow invariant manifold as shown in the previous section 
(compare to the method of slow invariant manifold for describing kinetics of dissipative 
systems |17j). The merits of the approach with different scales that simplifies both the 
procedure of finding the admitted group and construction of the group invariant solutions 
point to the quest for future potential applications. 
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